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Abstract 

We present calculations of the beauty photoproduction at HERA collider in the frame- 
work of the /cy-factorization approach. Both direct and resolved photon contributions are 
taken into account. The unintegrated gluon densities in a proton and in a photon are ob- 
tained from the full CCFM, from unified BFKL-DGLAP evolution equations as well as from 
the Kimber-Martin-Ryskin prescription. We investigate different production rates (both in- 
clusive and associated with hadronic jets) and compare our theoretical predictions with the 
recent experimental data taken by the HI and ZEUS collaborations. Special attention is put 
on the x° bs variable which is sensitive to the relative contributions to the beauty production 
cross section. 



1 Introduction 

The beauty production at high energies is a subject of intensive study from both theo- 
retical and experimental points of view. First measurements [1] of the 6-quark cross sections 
at HERA were significantly higher than QCD predictions calculated at next-to-leading or- 
der (NLO). Similar observations were made in hadron-hadron collisions at Tevatron [2] and 
also in photon-photon interactions at LEP2 [3]. In last case, the theoretical NLO QCD 
predictions are more than three standard deviations below the experimental data. At Teva- 
tron, recent analisys indicates that the overall description of the data can be improved [4] 
by adopting the non-perturbative fragmentation function of the 6-quark into the I?-meson: 
an appropriate treatment of the 6-quark fragmentation properties considerably reduces the 
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disagreement between measured beauty cross section and the corresponding NLO QCD cal- 
culations. Recently HI and ZEUS collaborations have reported important data [5-7] on 
the beauty photoproduction (both inclusive and associated with hadronic jets) in electron- 
proton collisions at HERA which refer to small values of the Bjorken scaling variable x. 
These data are in a reasonable agreement with NLO QCD predictions or somewhat higher. 
Some disagreement is observed [7] mainly at small decay muon and/or associated jet trans- 
verse momenta. But the large excess of the first measurements over NLO QCD, reported [1] 
by the HI collaboration, is not confirmed. In the present paper to analyze the HI and ZEUS 
data we will apply the so-called /c T -factorization [8, 9] (or semi-hard [10, 11]) approach of 
QCD since the beauty production at HERA is dominated by the photon-gluon or gluon- 
gluon fusion (direct and resolved photon contributions, respectively) and therefore sensitive 
to the gluon densities in a proton and in a photon at small values of x. 

The /c T -factorization approach is based on the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [12] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [13] gluon evolution which are 
valid at small x since here large logarithmic terms proportional to In 1/x are summed up to 
all orders of perturbation theory (in the leading logarithmic approximation). It is in contrast 
with the popular Dokshitzer-Gribov-Lipatov-Altarelli-Parizi (DGLAP) [14] strategy where 
only large logarithmic terms proportional to ln/i 2 are taken into account. The basic dynam- 
ical quantity of the /c T -factorization approach is the so-called unintegrated (k T -dependent) 
gluon distribution A(x, k^, /i 2 ) which determines the probability to find a gluon carrying the 
longitudinal momentum fraction x and the transverse momentum ky at the probing scale 
fi 2 . The unintegrated gluon distribution can be obtained from the analytical or numerical 
solution of the BFKL or CCFM evolution equations. Similar to DGLAP, to calculate the 
cross sections of any physical process the unintegrated gluon density A(x, k 2 -, /z 2 ) has to be 
convoluted [8-11] with the relevant partonic cross section a. But as the virtualities of the 
propagating gluons are no longer ordered, the partonic cross section has to be taken off mass 
shell (k^-dependent). It is in clear contrast with the DGLAP scheme (so-called collinear 
factorization). Since gluons in initial state are not on-shell and are characterized by vir- 
tual masses (proportional to their transverse momentum), it also assumes a modification of 
their polarization density matrix [10, 11]. In particular, the polarization vector of a gluon is 
no longer purely transversal, but acquires an admixture of longitudinal and time-like com- 
ponents. Other important properties of the /^-factorization formalism are the additional 
contribution to the cross sections due to the integration over the k 2 ^ region above fi 2 and the 
broadening of the transverse momentum distributions due to extra transverse momentum of 
the colliding partons. 

Some applications of the /^-factorization approach supplemented with the BFKL and 
CCFM evolution to the 6-quark production at high energies were discussed in [15-26]. It 
was shown [17-21] that the beauty cross section at Tevatron can be consistently described 
in the framework of this approach. However, a substantial discrepancy between theory and 
experiment is still found [22-25] for the 6-quark production in 77 collisions at LEP2, not 
being cured by the /cy-factorization 1 . At HERA, the inclusive beauty photoproduction has 
been investigated [18, 22, 26] and comparisons with the first HI measurements [1] have been 
done. It was concluded that the /c^-factorization approach provide a reasonable description 

1 Some discussions of this problem may be found in [22, 24]. 
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of the data within the large theoretical and experimental uncertainties. In [18, 26] the Monte- 
Carlo generator Cascade [27] has been used to predict the cross section of the 6-quark and 
dijet associated photoproduction. However, all calculations [18, 22, 26] deal with the total 
cross sections only. A number of important differential cross sections (such as transverse 
momentum and pseudo-rapidity distributions) has not been considered and comparisons 
with the recent HI and ZEUS measurements [5-7] have not been made. 

In the present paper we will study the beauty production at HERA in more detail. We 
investigate a number of different photoproduction rates (in particular, the transverse mo- 
mentum and pseudo-rapidity distributions of muons which originate from the semi-leptonic 
decays of 6-quarks) and make comparisons with the recent HI and ZEUS data [5-7]. Both 
direct (7(7 — > bb) and resolved photon contributions (gg — > bb) will be taken into account 2 . 
Our analysis covers also both inclusive and dijet associated 6-quark production. In last case 
special attention will be put on the x° bs variable since this quantity is sensitive to the relative 
contributions to the cross section from different production mechanisms. One of the purposes 
of this paper is to investigate the specific /^-factorization effects in the &-quark production 
at HERA. In the numerical analysis we test the unintegrated gluon distributions which are 
obtained from the full CCFM, unified BFKL-DGLAP [31] evolution equations and from the 
conventional parton densities (using the Kimber-Martin-Ryskin prescription [32]). We would 
like to note that this study is the continuation of our previous investigations [24, 30] where 
we have discussed, in particular, the charm production at LEP2 [24] and HERA [30]. 

The outline of our paper is following. In Section 2 we recall the basic formulas of the 
/cT-factorization approach with a brief review of calculation steps. In Section 3 we present 
the numerical results of our calculations and a dicussion. Finally, in Section 4, we give some 
conclusions. 

2 Basic formulas 

2.1 Kinematics 

We start from the gluon-gluon fusion subprocess. Let p e and p p be the four-momenta of 
the initial electron and proton, k\ and k 2 the four-momenta of the incoming off-shell gluons, 
and pb and pi the four-momenta of the produced beauty quarks. In our analysis below we 
will use the Sudakov decomposition, which has the following form: 



where k±T, k 2 x, P bT an< ^ Pit are the transverse four-momenta of the corresponding particles. 
It is important that k 2 T = — k\ T 7^ and k| T = —k\ T ^ 0. If we make replacement k\ — > p e 
and set x\ — 1 and kyr = 0, then we easily obtain more simpler formulas corresponding to 
photon-gluon fusion subprocess. In the ep center-of-mass frame we can write 



2 The 6-quark excitation processes bg — > bg are automatically included in the fey-factorization approach, 
as it was demonstrated in [28-30]. 



p b = a x p e + fixpp + p bT , Pi = a 2 p e + fcPp + PbT, 
ki = xip e + k 1T , k 2 = x 2 p p + k 2T , 



(1) 



p e = v^/2(l, 0, 0, 1), Pp = v^/2(l, 0, 0,-1), 



(2) 
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where s = (p e + p p ) 2 is the total energy of the process under consideration and we neglect 
the masses of the incoming particles. The Sudakov variables are expressed as follows: 

m bT ( v m- bT 

(3) 

Pi = exp(-y & ), fa = exp(-j/5), 
V s V s 

where rribT and mi T are the transverse masses of the produced quarks, and y b and m are 
their rapidities (in the ep center-of-mass frame). From the conservation laws we can easily 
obtain the following conditions: 

x 1 =ai + a 2 , x 2 = (3i+(3 2 , ki T + k 2T = PbT + Pit- ( 4 ) 

The variable x° bs is often used in the analysis of the data which contain the dijet samples. 
This variable, which is the fraction of the photon momentum contributing to the production 
of two hadronic jets with transverse energies E^ 1 and E^ 2 , experimentally is defined [6, 7] 

as 



_ E^e-^ + E^e-^ 



where yE e is the initial photon energy and rf cti are the pseudo-rapidities of these jets. The 
pseudo-rapidities rf oti are defined as rf eti = — lntan(^ ct ^/2), where 9 ieti are the polar angles 
of the jets with respect to the proton beam. Note that the selection of x° bs > 0.75 and 
a;° bs < 0.75 yields samples enriched in direct and resolved photon processes, respectively. 

2.2 Cross section for beauty photoproduction 

The main formulas for the total and differential beauty production cross sections were 
obtained in our previous papers [19, 24]. Here we recall some of them. In general case, the 
cross section a according to /c^-factorization theorem can be written as a convolution 

/dz 
— dk^C(x/z,k^, fi 2 )A(x, k|,/i 2 ), (6) 

where C(x, k^,/i 2 ) is the coefficient function corresponding to relevant partonic subprocess 
under consideration. So, the direct photon contribution to the differential cross section of 
7J9 — > bb + X process is given by 

da^\ 1P ^bb + X) _ f \M\ 2 ( ig *^bb) dfadfa 

dy h d V \ T 'J 16ir(x 2S y(l - a ± ) A[ 2 ' 2T ^ } 2T 2n 2vr ' 10 

where |-M| 2 (7g* — > bb) is the squared off-shell matrix element which depends on the trans- 
verse momentum k| T , <p 2 and 0& are the azimuthal angles of the initial virtual gluon and 
the produced quark, respectively. The formula for the resolved photon contribution can be 
obtained by the similar way. But one should keep in mind that convolution in (6) should 
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be made also with the unintegrated gluon distribution A^(x, k 2 -.,/! 2 ) in a photon. The final 
expression for the differential cross section has the form 

da (lcs \ 1P ^bb + X) _ f \M\ 2 (g*g* -> 66) 



/ 



dy b dpl T J 167r(a;ia;2S; ^ 

i 1 2 2\ At 1 2 2\ ti 2 t, 2 j ^4>1 ^02 d(f) b 

xAy{xi,k 1T , fi )A{x 2 ,\l 2T ,h )dk 1T dk 2T dy- b , 

Z7T Z7T Z7T 



where 0i is the azimuthal angle of the initial virtual gluon having fraction X\ of a ini- 
tial photon longitudinal momentum. It is important that squared off-shell matrix element 
\M\ 2 (g*g* — > 66) depends on the both transverse momenta k 2 T and k| T . The analytic ex- 
pressions for the |-M| 2 (7g* — > 66) and | | 2 (f?*(7* — > 66) have been evaluated in our previous 
papers [19, 24]. Note that if we average (7) and (8) over and k 2 ^ and take the limit 
k 2 T — > and k| T — > 0, then we obtain well-known formulas corresponding to the leading- 
order (LO) QCD calculations. 

The recent experimental data [5-7] taken by the HI and ZEUS collaborations refer to 
the 6-quark photoproduction in ep collisions, where electron is scattered at small angle and 
the mediating photon is almost real (Q 2 ~ 0). Therefore 7^ cross sections (7) and (8) needs 
to be weighted with the photon flux in the electron: 

da(ep ^bb + X) = J f y/e (y)dyda( 1P ^bb + X), (9) 

where y is a fraction of the initial electron energy taken by the photon in the laboratory 
frame, and we use the Weizacker- Williams approximation for the bremsstrahlung photon 
distribution from an electron: 

Here a em is Sommerfeld's fine structure constant, m e is the electron mass, <5mm = m ly 2 / (1 ~~ 
y) 2 and Q 2 mx = lGeV 2 , which is a typical value for the recent photoproduction measure- 
ments at HERA. 

The multidimensional integration in (7), (8) and (9) has been performed by means of 
the Monte Carlo technique, using the routine Vegas [33]. The full C++ code is available 
from the authors on request 3 . This code is practically identical to that used in [30], with 
exception that now we apply it to calculate beauty production instead charm. 



3 Numerical results 

We now are in a position to present our numerical results. First we describe our theoret- 
ical input and the kinematical conditions. 
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3.1 Theoretical uncertainties 



There are several parameters which determined the normalization factor of the cross 
sections (7) and (8): the beauty mass m&, the factorization and normalisation scales lif 
and [xr and the unintegrated gluon distributions in a proton A(x, k^, fx 2 ) and in a photon 
^l 7 (a;,k^/i 2 ). 

Concerning the unintegrated gluon densities in a proton, in the numerical calculations 
we used five different sets of them, namely the J2003 (set 1 - - 3) [21], KMS [31] and 
KMR [32]. All these distributions are widely discussed in the literature (see, for example, 
review [34, 35] for more information). Here we only shortly discuss their characteristic 
properties. First, three sets of the J2003 gluon density have been obtained [21] from the 
numerical solution of the full CCFM equation. The input parameters were fitted to describe 
the proton structure function F 2 (x, Q 2 ). Note that the J2003 set 1 and J2003 set 3 densities 
contain only singular terms in the CCFM splitting function P gg (z). The J2003 set 2 gluon 
density takes into account the additional non-singlular terms 4 . These distributions have 
been applied in the analysis of the forward jet production at HERA and charm and bottom 
production at Tevatron [21] (in the framework of Monte-Carlo generator Cascade) and 
have been used also in our calculations [30]. 

Another set (the KMS) [31] was obtained from a unified BFKL-DGLAP description 
of F 2 (x, Q 2 ) data and includes the so-called consistency constraint [36]. The consistency 
constraint introduces a large correction to the LO BFKL equation. It was argued [37] that 
about 70% of the full NLO corrections to the BFKL exponent A are effectively included in 
this constraint. The KMS gluon density is successful in description of the beauty production 
at Tevatron [17, 19] and J/ip meson photo- and leptoproduction at HERA [38, 39]. 

The last, fifth unintegrated gluon distribution A(x, k^, fx 2 ) used here (the so-called KMR 
distribution) is the one which was originally proposed in [32]. The KMR approach is the 
formalism to construct unintegrated gluon distribution from the known conventional parton 
(quark and gluon) densities. It accounts for the angular-ordering (which comes from the 
coherence effects in gluon emission) as well as the main part of the collinear higher-order 
QCD corrections. The key observation here is that the fx dependence of the unintegrated 
parton distribution enters at the last step of the evolution, and therefore single scale evolution 
equations (DGLAP or unified BFKL-DGLAP) can be used up to this step. Also it was 
shown [32] that the unintegrated distributions obtained via unified BFKL-DGLAP evolution 
are rather similar to those based on the pure DGLAP equations. It is because the imposition 
of the angular ordering constraint is more important [32] than including the BFKL effects. 
Based on this point, in our further calculations we use much more simpler DGLAP equation 
up to the last evolution step 5 . Note that the KMR parton densities in a proton were used, 
in particular, to describe the prompt photon photo- and hadroproduction at HERA [41] and 
Tevatron [42, 43]. 

In the case of a real photon, we have tested two different sets of the unintegrated gluon 
densities Aj(x, k^, li 2 ). First of them was obtained [23] from the numerical solution of the 
full CCFM equation (which has been also formulated for the photon). Here we will use this 
gluon density together with the three sets of the J2003 distribution when calculating the 

4 See Ref. [21] for more details. 

5 We have used the standard GRV (LO) parametrizations [40] of the collinear quark and gluon densities. 
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resolved photon contribution (8). Also in order to obtain the unintegrated gluon density in 
a photon we will apply the KMR method to the standard GRV parton distributions [40]. In 
the numerical calculations we will use it together with the KMR distributions in a proton. 
Note that both gluon densities *4. 7 (x, k^, fx 2 ) discussed here have been already applied in the 
analysis of the heavy (charm and beauty) quark [22-25] and J/ip meson [24, 25] production 
in 77 collisions at LEP2. 

We would like to point out that at present there is not the unintegrated gluon distribution 
corresponding to the unified BFKL-DGLAP evolution in a photon. Therefore we will not 
take into account the resolved photon contribution (8) in the case of KMS gluon. 

Also the significant theoretical uncertainties in our results connect with the choice of the 
factorization and renormalization scales. First of them is related to the evolution of the gluon 
distributions, the other is responsible for the strong coupling constant a s (fi 2 R ). The optimal 
values of these scales are such that the contribution of higher orders in the perturbative 
expansion is minimal. As it often done [11, 15-26, 44] for beauty production, we choose 

the renormalization and factorization scales to be equal: lir = Lip = li = \Jm 2 + (p^), 
where (pf.) is set to the average pf, of the beauty quark and antiquark 6 . Note that in the 
present paper we concentrate mostly on the non-collinear gluon evolution in the proton and 
do not study the scale dependence of our results. To completeness, we take the 6-quark mass 
m b = 4.75 GeV and use LO formula for the coupling constant a s (/i 2 ) with rif = 4 active 
quark flavours at A QCD = 200 MeV, such that cc s (Mf) = 0.1232. 

3.2 Inclusive beauty photoproduction 

The recent experimental data [5-7] for the inclusive beauty photoproduction at HERA 
comes from both the HI and ZEUS collaborations. The 6-quark total cross section for 
Pt > Pt" 1 as wen as the beauty transverse momentum distribution have been determined. 
The ZEUS data [5, 6] refer to the kinematical region 7 defined by \rj b \ < 2 and Q 2 < 1 GeV 2 , 
where rj b is the beauty pseudo-rapidity. The fraction y of the electron energy transferred to 
the photon is restricted to the range 0.2 < y < 0.8. 

In Figs. 1 and 2 we show our predictions in comparison to the ZEUS data [5, 6]. The 
solid, dashed, dash-dotted, dotted and short dash-dotted curves correspond to the results 
obtained with the J2003 set 1 — 3, KMR and KMS unintegrated gluon densities, respec- 
tively. One can see that overall agreement between our predictions and experimental data is 
a very good. All three sets of the J2003 distribution as well as the KMS gluon density give 
results which are rather close to each other (except large p b T region where the KMS density 
predicts more hard behaviour). We find also a some enhancement of the estimated cross 
sections as compared with the collinear NLO QCD calculations which lie somewhat below 
the measurements but still agree with the data within the scale uncertainties. This enhance- 
ment comes, in particular, from the non-zero transverse momentum of the incoming off-shell 
gluons. Note that the KMR gluon distribution gives results which lie below the ZEUS data 
and which are very similar to the NLO QCD predictions. This observation coincides with 

6 We use special choice — in the case of KMS gluon, as it was originally proposed in [31]. 
7 Here and in the following all kinematic quantities are given in the laboratory frame where positive OZ 
axis direction is given by the proton beam. 
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Source 


a(ep — > e'bb + X) [nb] 


HI measurement [1] 
Cascade [26] 

Torino c*+ 1 
jzuuo sei i 

J2003 set 2 

J2003 set 3 

KMR 

KMS 


14.8 ± 1.3 (stat.) (sys.) 

r 9 +l.l 
U. ( o 

6.62 
7.16 
3.91 
7.57 



Table 1: The total cross section of the inclusive beauty photoproduction in electron-proton 
collisions at Q 2 < 1 GeV 2 . 

the ones [41]. Such underestimation can be explained by the fact that leading logarithmic 
terms proprtional to In 1/x are not included into the KMR formalism. 

Also the total inclusive beauty cross section a(ep — > ebb + X) has been measured [1] by 
the HI collaboration and it was found to be equal to 14.8 ± 1.3 (stat.) j^J (sys.) nb for 
Q 2 < 1 GeV 2 . The collinear NLO QCD calculations predict a cross section which is about a 
factor of 4 below the HI measurements [1] . The results of our calculations supplemented with 
the different unintegrated gluon densities are collected in Table 1. Also the predictions of the 
Monte-Carlo generator Cascade [26] are shown for comparison. One can see that earlier 
HI data [1] exceed our theoretical estimations by a factor about 2. However, recent analysis 
which has been performed in [6, 7] does not confirm the results of the first measurements [1]. 
So, the cross section for muon coming from b decays in dijet photoproduction events was 
found [7] to be significantly lower than one reported in [1]. Therefore we can expect that the 
inclusive 6-quark cross section (which can be obtained after extrapolation of dijet and muon 
cross section to the full phase space) will be reduced and agreement with our predictions 
will be significantly improved. 

In general, we can conclude that the cross sections of inclusive beauty photoproduc- 
tion calculated in the /c^-factorization formalism (supplemented with the CCFM or unified 
BFKL-DGLAP evolution) are larger by 30 — 40% than ones calculated at NLO level of 
collinear QCD. Our results for the total and differential cross sections are in a better agree- 
ment with the HI and ZEUS data than the NLO QCD predictions. We find also that the 
individual contributions from the photon-gluon and gluon-gluon fusion to the inclusive b- 
quark cross section in the /cy-factorization approach is about 85 and 15%, respectively. This 
is in agreement with the results presented in [22] where the KMR and GBW unintegrated 
gluon densities has been used. 

3.3 Dijet associated beauty photoproduction 

Now we demonstrate how the /^-factorization approach can be used to calculate the 
semi-inclusive beauty photoproduction rates. The basic photon-gluon or gluon-gluon fusion 
subprocesses give rise to two high-energy 6-quarks, which can further evolve into hadron 
jets. In our calculations the produced quarks (with their known kinematical parameters) 
were taken to play the role of the final jets. These two quarks are accompanied by a number 
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of gluons radiated in the course of the gluon evolution. As it has been noted in [28], on the 
average the gluon transverse momentum decreases from the hard interaction block towards 
the proton. We assume that the gluon emitted in the last evolution step and having the 
four-momenta k' compensates the whole transverse momentum of the gluon participating in 
the hard subprocess, i.e. k' T ~ — k T . All the other emitted gluons are collected together 
in the proton remnant, which is assumed 8 to carry only a negligible transverse momentum 
compared to k' T - This gluon gives rise to a final hadron jet with Ej. = |k' T | in addition to 
the jet produced in the hard subprocess. From these three hadron jets we choose the two 
ones carrying the largest transverse energies, and then compute the beauty and associated 
dijet photoproduction rates. 

The recent experimental data [6, 7] on the beauty and associated dijet production at 
HERA come from both HI and ZEUS collaborations. The ZEUS data [6] refer to the 
kinematical region defined by 0.2 < y < 0.8, Q 2 < 1 GeV 2 and given for jets with 1 > 
7 GeV, pp > 6 GeV and | r/ et | < 2.5. The measured cross sections have been presented for 
muons coming from semileptonic b decays in dijet events with p%< > 2.5 GeV and —1.6 < 
rj^ < 2.3. The more recent HI data [7] refer to the same kinematical region except another 
muon pseudo-rapidity requirement: —0.55 < rf < 1.1. To produce muons from 6-quarks 
in our theoretical calculations, we first convert 6-quarks into 5-hadrons using the Peterson 
fragmentation function [45] and then simulate their semileptonic decay according to the 
standard electroweak theory 9 . Our default set of the fragmentation parameter is e& = 0.0035. 

So, the transverse momentum and pseudo-rapidity distributions of the 6-quark decay 
muon for different kinematical region are shown in Figs. 3 — 6 in comparison to the HERA 
data. One can see that calculated cross sections (using the J2003 and KMS unintegrated 
gluon densities) agree very well with the experimental data except the low pt^ region (p^ < 
3 GeV) in Fig. 4. Note, however, that the behaviour of measured cross sections in this 
region is very different from each other in the HI and ZEUS data. The pj* distribution 
measured [7] by the HI collaboration falls steeply with increasing transverse momentum 
Pj,. The similar situation is observed also for the cross section measured as a function of 
the transverse momentum of leading jet (see Fig. 7). Our calculations give a less steep 
behaviour and are lower than the HI data in the lowest momentum bin by a factor of 2.5. 
But at higher transverse momenta p^ better agreement is obtained. Note that in the case of 
Py* distribution the discrepancy at low p^ is smaller, is about 1.5 times only. In contrast, 
a good description of the ZEUS data [6] (both in normalization and shape) for all values of 
Pj. is observed (see Fig. 3). Therefore there is some inconsistency between the data. 

Also the ZEUS collaboration have presented the data on the transverse momentum and 
pseudo-rapidity distributions of the jets associated with the muon (so-called //-jet) or B- 
hadron (6-jet). These jets reproduce the kinematic of the b (or b) quark in a good ap- 
proximation. The //-jet is defined as the jet containing the 5-hadron that decays into the 
muon. Similarly, the fo-jet is defined as the jet containing the B (or B) hadron. In Figs. 8 
- 11 we show our predictions for the transverse momentum and pseudo-rapidity distribu- 
tions of the //-jet and 6-jet in comparison to the ZEUS measurements [6]. One can see that 

8 Note that such assumption is also used in the KMR formalism. 

9 Of course, the muon transverse momenta spectra are sensitive to the fragmentation functions. However, 
this dependence is expected to be small as compared with the uncertainties coming from the unintegrated 
gluon densities in a proton and in a photon. 
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J2003 and KMS gluon densities give results which agree well with the data, although slightly 
overestimate the data at low j9j.~ jet (see Fig. 8). 

We would like to note that the KMS gluon provides a more hard transverse momentum 
distribution of the final muon (or jets) as compared with other unintegrated densities under 
consideration. Similar effect we have observed in the case of the inclusive beauty photopro- 
duction (in a previous section). Another interesting observation is that the dotted curves 
which obtained using the KMR unintegrated gluon lie below the HI and ZEUS data every- 
where. This fact confirms the assumption which was made in [41] that the KMR formalism 
results in some underestimation of the calculated cross sections. Also it is interesting that 
the difference in normalization between the KMS and J2003 predictions is rather small, is 
about 20% only. However, it is in the contrast with the D* and dijet associated photoproduc- 
tion at HERA which has been investigated in our previous paper [30], where we have found 
a relative large enhancement of the cross sections calculated using the KMS gluon density. 
The possible explanation of this fact is that the large 6-quark or J ftp meson mass (which 
provide a hard scale) make predictions of the perturbation theory of QCD more applicable. 

Next we concentrate on the very interesting subject of study which is connected with 
the individual contributions from the direct and resolved photon mechanisms to the cross 
section in the /c^-factorization approach. As it was already mentioned above, the i° bs variable 
(which corresponds at leading order to the fraction of the exchanged photon momentum in 
the hard scattering process) provides a tool to investigate the relative importance of different 
contributions. In LO approximation, direct photon events at parton level have x° bs ~ 1, 
while the resolved photon events populate the low values of x° bs . The same situation is 
observed in a NLO calculations, because in the three parton final state any of these partons 
are allowed to take any kinematically accessible value. In the A; T -factorization formalism the 
hardest transverse momentum parton emission can be anywhere in the evolution chain, and 
does not need to be closest to the photon as required by the strong /z 2 ordering in DGLAP. 
Thus, if the two hardest jets are produced by the bb pair, then x° bs is close to unity, but if a 
gluon from the initial cascade and one of the final 6-quarks form the two hardest transverse 
momentum jets, then x° bs < I. This statement is clearly demonstrated in Fig. 12 where 
separately shown the contributions from the photon-gluon (dashed curve) and gluon-gluon 
fusion (dash-dotted curve) subprocesses. The solid curve represents the sum of both these 
contributions. We have used here the KMR unintegrated gluon density for illustration. As 
it was expected, the gluon-gluon fusion events (with a gluon coming from the photon) are 
distributed over the whole x° hs range. It is clear that these events play important role at 
small values of a;° bs . Next, in agreement with the expectation for direct photon processes, the 
peak at high values of the a;° bs is observed. However, one can see that off-shell photon-gluon 
fusion results also in substantial tail at small values of x° bs . The existence of this plateau 
in the collinear approximation of QCD usually is attributted to the heavy quark excitation 
from resolved photon. In the /cT-factorization approach such plateau indicates the fact that 
the gluon radiated from evolution cascade appears to be harder than 6-quarks (produced in 
hard parton interaction) in a significant fraction of events [28-30]. 

In Fig. 13 and 14 we confront the x° bs distributions calculated in different kinematical 
regions with the HERA data [6, 7]. One can see that the J2003 and KMR unintegrated 
gluon densities give a reasonable description of the data but tend to slightly underestimate 
them at middle and low x° hs . In the case of KMS gluon this discrepancy is more significant 
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Source 


o(ep — > ebb + X — > ejj/i + X') [pb] 


HI measurement [7] 
NLO QCD (FMNR) [44] 
Cascade (J2003 set 2) 
Pythia [46] 
J2003 set 1 

Torino apt 9 
jzuuo sei z 

J2003 set 3 

KMR 

KMS 

KMS (m b = 4.5 GeV, A QCD = 250 MeV) 


38.4 ± 3.4 (stat.) ± 5.4 (sys.) 
23.81^ 
22.6 
20.9 
28.37 

29.25 
17.43 
33.87 
38.84 



Table 2: The total cross section of the beauty and associated dijet photoproduction obtained 
in the kinematic range -0.55 < jf < 1.1, p£ > 2.5 GeV, Q 2 < 1 GeV 2 , 0.2 < y < 0.8, 
p^ 1 > 7 GeV, p J ^ 2 > 6 GeV and \^ ct \ < 2.5. 

since the contribution from the gluon-gluon fusion subprocess are not taken into account 
here. The Monte-Carlo generator Cascade [27] (which generates low x° bs events via initial 
state gluon radiation without using a gluon density in a photon) also underestimate [7] the 
cross sections at low :r° bs . Note that the shapes of x° bs distributions predicted by the J2003 
and KMR densities differs from each other. This fact is connected with different properties 
of corresponding unintegrated gluon distributions in a proton and in a photon. In general, 
from Fig. 12 — 14 we can conclude that the gluon-gluon fusion contribution is important in 
description of the experimental data and that the behaviour of calculated x° bs distributions 
at low values of x° bs is strongly depends on the unintegrated gluon densities used. However, 
our calculations still reasonable agree with the HI and ZEUS data within the theoretical 
and experimental uncertainties. 

Now we turn to the total cross section of 6-quark and associated dijet photoproduction. 
The ZEUS collaboration has presented [6] the results for forward, barrel and rear muon- 
chambers regions which defined by —1.6 < rf < —0.9, p%< > 2.5 GeV (rear), —0.9 < i]^ < 1.3, 
p^ > 2.5 GeV (barrel) and 1.48 < rf < 2.3, p^ T > 2.5 GeV, p» > 4 GeV (forward). In Fig. 15 
we display the results of our calculations in comparison to the recent ZEUS data. One can 
see that our predictions (supplemented with the J2003 and KMS gluon densities) agree well 
with the measurements in the barrel region but underestimate the data in rear and forward 
ones. The main discrepancy is found in forward kinematical region where our predictions 
are below the data by a factor of 1.5. In Table 2 we compare the calculated cross section 
with the HI data [7] which has been obtained in another kinematical region (defined above). 
The predictions of Monte-Carlo programs Pythia [46], Cascade [27] as well as NLO QCD 
calculations (FMNR) [44] are also shown for comparison. Note that these results are in a good 
agreement with each other but are about 1.5 standard deviations below [7] the data. Our 
predictions are somewhat higher but still lie below the data, too. However, this discrepancy 
is not dramatic, because some reasonable variations in beauty mass m^, energy scale ji 2 or 
Aqcd parameter, namely 4.5 < m b < 5 GeV, /ig/2 < /i 2 < 2/Xq (where /i is the transverse 
mass of produced 6-quark) and 150 < Aq C d < 250 MeV, can completely eliminate the visible 
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disagreement. To be precise, we have repeated our calculations using the KMS gluon density 
with the m b = 4.5 GeV and Aq CD = 250 MeV. We obtained the value a = 38.84 pb which 
is close to the experimental data point a = 38.4 ± 3.4 ± 5.4 pb. 

Finally, we would like to note that in according to the analysis [6, 7] which was done by 
the HI and ZEUS collaborations, in order to obtain a realistic comparison of the data and 
NLO QCD calculations the corrections for hadronisation should be taken into account in the 
predictions. The correction factors are typically 0.8 — 1.1 depending on a bin [6, 7]. These 
factors are not accounted for in our analysis. 

4 Conclusions 

We have investigated the beauty production in electron-proton collisions at HERA in the 
fc^-factorization QCD approach. The different photoproduction rates (both inclusive and 
associated with hadronic jets) have been studied. We took into account both the direct and 
resolved photon contribution. In numerical analysis we have used the unintegrated gluon 
densities which are obtained from the full CCFM, from unified BFKL-DGLAP evolution 
equations (KMS) as well as from the Kimber-Martin-Ryskin prescription. Our investiga- 
tions were based on LO off-mass shell matrix elements for the photon-gluon and gluon-gluon 
fusion subprocesses. Special attention has been drawn to the x° bs variable since this quantity 
is sensitive to relative contributions to the cross section from the different production mech- 
anisms. We demonstrate the importance of gluon-gluon fusion subprocess in description of 
the experimental data at low values of x° bs . 

We have shown that the /^-factorization approach supplemented with the CCFM or 
BFKL-DGLAP evolved unintegrated gluon distributions (the J2003 or KMS densities) re- 
produces well the numerous HERA data on beauty production. At the same time we have 
obtained that the Kimber-Martin-Ryskin formalism results in some underestimation of the 
cross sections. This shows the importance of a detail understanding of the non-collinear 
parton evolution process. 
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Figure 1: The inclusive beauty cross section as a function of p™ m calculated at |t7 & | < 2, Q 2 < 
1 GeV 2 and 0.2 < y < 0.8. The solid, dashed, dash-dotted, dotted and short dash-dotted 
curves correspond to the J2003 set 1 — 3, KMR and KMS unintegrated gluon distributions, 
respectively. The experimental data are from ZEUS [5]. 
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Figure 2: The differential cross section da/dp b T for the inclusive beauty production calculated 
at \r] h \ < 2, Q 2 < 1 GeV 2 and 0.2 < y < 0.8. All curves are the same as in Fig. 1. The 
experimental data are from ZEUS [5, 6]. 



16 




4 6 8 10 

Pf (GeV) 

Figure 3: The differential cross section da/dp^ for dijets with an associated muon coming 
from b decays in the kinematic range —1.6 < T]^ < 2.3, Q 2 < 1 GeV 2 , 0.2 < y < 0.8, 
p 1 ? 1 > 7 GeV, p^ 2 > 6 GeV and |r? jct | < 2.5. All curves are the same as in Fig. 1. The 
experimental data are from ZEUS [6] . 
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Figure 4: The differential cross section da/dp^ for dijets with an associated muon coming 
from b decays in the kinematic range —0.55 < < 1.1, Q 2 < 1 GeV 2 , 0.2 < y < 0.8, 
p^ 1 > 7 GeV, p S T 2 > 6 GeV and |r] jct | < 2.5. All curves are the same as in Fig. 1. The 
experimental data are from HI [7]. 
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Figure 5: The differential cross section da jdrf for dijets with an associated muon coming 
from b decays in the kinematic range p% > 2.5 GeV, Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p 3 ^ 1 > 
7 GeV, p S T 2 > 6 GeV and |?? jct | < 2.5. All curves are the same as in Fig. 1. The experimental 
data are from ZEUS [6]. 
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Figure 6: The differential cross section do jdrf for dijets with an associated muon coming 
from b decays in the kinematic range pj* > 2.5 GeV, Q 2 < 1 GeV 2 , 0.2 < y < 0.8, 1 > 
7 GeV, p S T 2 > 6 GeV and |?? jct | < 2.5. All curves are the same as in Fig. 1. The experimental 
data are from HI [7]. 
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Figure 7: The leading jet transverse momentum distribution da/dpj? for dijets with an 
associated muon coming from b decays in the kinematic range —0.55 < rf < 1.1, Q 2 < 
1 GeV 2 , 0.2 < y < 0.8, p^ 1 > 7 GeV, p^ 2 > 6 GeV and |^ ct | < 2.5. All curves here are the 
same as in Fig. 1. The experimental data are from HI [7]. 
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Figure 8: The transverse momentum distribution of the jet associated to the muon coming 
from b decays in the kinematic range —0.55 < t/ 1 < 1.1, > 2.5 GeV, Q 2 < 1 GeV 2 , 
0.2 < y < 0.8 and |^^J ct | < 2.5. All curves are the same as in Fig. 1. The experimental 
data are from ZEUS [6]. 
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Figure 9: The pseudo-rapidity distribution of the jet associated to the muon coming from b 
decays in the kinematic range -0.55 < -rf < 1.1, p%> > 2.5 GeV, Q 2 < 1 GeV 2 , 0.2 < y < 0.8 
and p£" jet > 6 GeV. All curves are the same as in Fig. 1. The experimental data are from 
ZEUS [6]. 
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Figure 10: The transverse momentum distribution of the jet containing a 5-hadron in the 
kinematic range -0.55 < rf < 1.1, > 2.5 GeV, Q 2 < 1 GeV 2 , 0.2 < y < 0.8 and 
1^6— jet | < 2_5_ All curves are the same as in Fig. 1. The experimental data are from 
ZEUS [6]. 
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Figure 11: The pseudo-rapidity distribution of the jet containing a 5-hadron in the kinematic 
range -0.55 < rf < 1.1, Vt > 2 - 5 GeV > Q 2 < 1 GeV2 > 0.2 < y < 0.8 and |^~j ct | < 2.5. All 
curves are the same as in Fig. 1. The experimental data are from ZEUS [6]. 
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Figure 12: The differential cross section da/dx^ 38 for dijets with an associated muon coming 
from b decays in the kinematic range —1.6 < nf 1 < 2.3, pj. > 2.5 GeV, Q 2 < 1 GeV 2 , 
0.2 < y < 0.8, p^ 1 > 7 GeV, pfi 2 > 6 GeV and |^ ct | < 2.5. Separately shown the 
contributions from the photon-gluon (dashed curve) and gluon-gluon fusion (dash-dotted 
curve). Solid curve represents the sum of both these contributions. The KMR unintegrated 
gluon densities in a proton and in a photon has been used. 
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Figure 13: The differential cross section da/dx^ for dijets with an associated muon coming 
from b decays in the kinematic range —1.6 < t/ 4 < 2.3, > 2.5 GeV, Q 2 < 1 GeV 2 , 
0.2 < y < 0.8, pf* 1 > 7 GeV, pf 2 > 6 GeV and |r? jct | < 2.5. All curves are the same as in 
Fig. 1. The experimental data are from ZEUS [6]. 
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Figure 14: The differential cross section da/dx^ 38 for dijets with an associated muon coming 
from b decays in the kinematic range —0.55 < rf < 1.1, > 2.5 GeV, Q 2 < 1 GeV 2 , 
0.2 < y < 0.8, pf 1 > 7 GeV, p^ 2 > 6 GeV and |?? jct | < 2.5. All curves are the same as in 
Fig. 1. The experimental data are from HI [7]. 
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Figure 15: The cross section for muon coming from b decays in dijet events calculated in the 
rear, barrel and forward kinematical regions (see text). The cuts are applied: Q 2 < 1 GeV 2 , 
0.2 < y < 0.8, p S T h > 7 GeV, pi^ 2 > 6 GeV and |^ jct | < 2.5. All curves are the same as in 
Fig. I. The experimental data are from ZEUS [6]. 
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